International Journal of Public Health Science (IJPHS) 
Vol. 10, No. 2, June 2021, pp. 231~240 
ISSN: 2252-8806, DOI: 10.11591/ijphs.v10i12.20713 Oj pe | 


A spatiotemporal analysis of dengue hemorrhagic fever in 
Banyumas, Indonesia 


Dwi Sarwani Sri Rejeki’, Nunung Nurhayati’, Budi Aji° 
'’Department of Public Health, Faculty of Health Sciences, Universitas Jenderal Soedirman, Indonesia 
"Department of Mathematics, Faculty of Mathematics and Science, Universitas Jenderal Soedirman, Indonesia 


Article Info ABSTRACT 

Article history: Dengue hemorrhagic fever (DHF) is a public health problem in the world and 
also in Indonesia. One of the districts in Central Java that is still having 
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dengue fever patients in 2018 using as many as 57 cases. Furthermore, the 
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Risk zones data analysis used includes overlay, buffering, and clustering with SaTScan 
Spatiotemporal and ArcGis software. The results showed that there was a clustering of 
Transmission dengue cases in Banyumas, with one primary and three secondary clusters 
Variation detected. The primary cluster occurred in March-April 2018, involving four 


sub-districts in urban areas. It was then concluded that the significantly 
identified clusters indicate a transmission of dengue fever in the Banyumas 
area with a radius of three kilometers. 
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1, INTRODUCTION 

Dengue hemorrhagic fever (DHF) is a vector-borne disease that is widespread in the world, 
especially in tropical and subtropical climates. This infection is caused by the dengue virus which is 
transmitted from Aedes aegpty or Aedes albopictus mosquito species. According to WHO, Asia Pacific bore 
75% of the world's dengue burden between 2004 and 2010, while Indonesia was reported as the 2nd country 
with the largest cases among 30 endemic countries [1]. The Indonesian Ministry of Health recorded dengue 
fever incidence rate (IR) of 50.75 per 100,000 populations in 2015, with a case fatality rate (CFR) of 0.82%. 
This figure increased in 2016 with an IR of 79.92 per 100,000 populations and a CFR of 0.78%. However, 
the number of cases decreased in 2017, with IR recorded at 26.12 per 100,000 population and a CFR of 
0.72% [2]. 

Banyumas is a contributor to dengue cases in Central Java. The number of cases in Banyumas was 
264 in 2015 (IR=16.43), 990 in 2016 (IR=55.07), and 68 in 2017 (IR=3.35). In 2018, DHF cases decreased to 
55 (IR 2.75 per 100,000) and the CFR value was 3.64% [3]. The results of a qualitative study of DHF control 
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program in Banyumas showed that the input components of the prevention and control program (funds, 
methods, and time) are relatively sufficient. However, the implementation of programs such as fogging focus, 
larvicide, regular eradication of mosquito nests is still not optimal [4]. 

Geospatial analysis is useful in controlling vector-borne diseases such as dengue and malaria. This 
analysis is a spatial aspect that shows the location and position of an object or event that is below, on, or 
above the earth's surface which is stated in a certain coordinate system [5]. Changing spatial data into 
geospatial information requires a geographical information system (GIS) and spatial analysis. Furthermore, 
the Geographic information systems are to determine the analysis of local variations and geographic risk 
zones, as well as to measure disease control interventions [6]. The combination of spatial and time aspects is 
called spatiotemporal. The results of the systematic review of dengue fever showed a useful spatial analysis 
in the public health sector, including identifying the risk of the fever. This method is used to advance 
scientific knowledge about the dynamics of dengue transmission, surveillance strategies, and to reduce the 
social costs associated with the transmission [7]. 

The spatial analysis of DHF in Putrajaya, Malaysia showed that the incidence of dengue occurs in 
clusters with an average radius of 264 meters [8]. Furthermore, findings in Taiwan showed that there is a 
spatial grouping of DHF cases in metropolitan areas and Southern Taiwan which is a green space for 
landscaping [9]. In Queensland, Australia, the incidence of dengue fever also clustered around Cairns and 
occurred in autumn [10]. The spatial and temporal results of the fever in Guangdong China showed the 
incidence of centralized outbreaks in the Yuexiu District and spread to several surrounding districts in 
Guangzhou and Guangdong Province [11]. Provinces with a large number of case counts, urban population, 
urban land cover and incoming flight passengers are associated to the epidemic prone cluster of dengue in 
Thailand [12]. From several spatial studies, it is important to pinpoint the area of distribution and timing in 
controlling DHF. 

GIS and spatial statistics were frequently used to characterize spatiotemporal patterns of dengue and 
other infectious diseases. Acharya et al. in Nepal assessed spatiotemporal patterns of dengue in Nepal 
combining the address and district levels [13]. SaTScan software provides a powerful tool to detect, 
delineate, and validate disease clusters, risk population, and factors associated with them over space and 
time. Further SaTScan adjusts for confounding variables, and reduces preselection bias regarding the size and 
location of cluster [13]. This study uses analysis spasiotemporal at the subdistrict level. Spatiotemporal 
analysis indices correlated with the highest DHF density: high transmission intensity only, long duration, and 
high transmission intensity plus long duration [14]. This analysis could provide useful information to identify 
areas at high risk of DHF and time. This information about analysis spasiotemporal DHF in subdistrict level 
may be helpful to the development of related epidemic prevention and control strategies in Banyumas 
district. 

The results of the spatiotemporal information on DHF are used as a tool in formulating policies, 
making decisions, and implementing activities related to DHF. Also, spatiotemporal analysis in epidemiology 
is very useful, especially to evaluate the occurrence of differences in incidence by geographic area and 
identify disease clustering. The benefits of clustering analysis include showing the geographical surveillance 
of disease, identifying disease clusters spatially or space-time and, determining whether the clusters are 
statistically significant [15]. The results of this spatiotemporal analysis are used as control information, 
therefore the target and timing of its implementation are right. Therefore, this study aimed to determine the 
spatiotemporal distribution and grouping of DHF cases in sub district level in Banyumas 2018. 


2. RESEARCH METHOD 

The study design was observational with a cross-sectional spatial analysis. It was conducted in the 
Banyumas area with 57 DHF cases recorded at the Banyumas Health Office in 2018. The studied variables 
include the house of dengue fever patients, date of illness, location of health center services, population 
density in villages or sub-districts, land use, rain, rivers, and rock types. The data were collected from several 
sources, namely primary and secondary. The secondary data were obtained from the public health center 
(Puskesmas), Statistics Indonesia, National Meteorology Climatology Geophysics Agency, Civil Registry 
Service Office, and Regional Planning Organization. Furthermore, the primary data which include the 
coordinates of the patients’ houses, the coordinates of Puskesmas, and height of the house were obtained 
using the geographical positioning system (GPS) tool. The geographic coordinates of the case are the 
coordinates of the house, which was done outside and in front of the house. The analysis includes overlay 
analysis, buffering, and clustering. Overlay analysis was used to determine the spread of DHF cases based on 
several variables such as land use, population density, rainfall, and rock types. Meanwhile, population 
density in this study is the number of people in each village. This variable was obtained by dividing the total 
population by the area of each village (km*). Furthermore, buffering analysis was used to determine the 
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spread of the disease, and the spatial elements used in this analysis are the existence of health services or 
centers and rivers. Clustering analysis was conducted with SaTScan software to determine the cluster or 
grouping of dengue events in space and time (space-time permutation model). Analysis of the space-time 
permutation model was visualized using the Arcmap 10.1 program. In addition, SaTScan analysis used the 
Space-Time Permutation Model (Likelihood Ratio Test) by analyzing significant clusters, namely the most 
likely and secondary cluster. 


3. RESULTS 

Banyumas Regency is a lowland with an average height of +108 meters above sea level, and located 
between 70 15 '05 "- 70 37' 10" south latitude and 1080 39 '17 "—1090 27" 15 " north latitude. This regency 
has a land area of 1,327.59 km’ consisting of 27 sub-districts, with Cilongok as the largest (105.34 km’) and 
West Purwokerto as the smallest sub-district (7.40 km’). Based on the elevation (height above sea level), the 
plains in Banyumas Regency consist of: Om-100m=54.86%, 101m-500m=45.14%. This regency has a 
tropical climate with two seasons, namely dry and rainy. Throughout 2017, there was rainfall which 
fluctuated for 153 days and varied by month. The highest rainfall was recorded in October at 648.73mm, 
while the lowest was in August at 12.45 mm. Furthermore, Banyumas Regency has a wet tropical climate 
with an average temperature of 26.5°C. The minimum temperature is around 24.4°C and the maximum is 
around 30.9°C [16]. 

Banyumas has 27 sub-districts with 331 villages or wards and is led by a Regent and Deputy 
Regent. Based on population projections in 2017, the number of people in this regency is 1,665,025, 
consisting of 831,816 men and 833,209 women. Cilongok is the sub-district with the largest population with 
118,366 people and Purwojati is the smallest with 32,025 people. The population density in Banyumas 
Regency in 2017 was 1,254 people/km’ with an average of 3.6 people per household. The population density 
in 27 sub-districts is quite diverse, with the highest in North Purwokerto of 7,326 people/km’ and the lowest 
in Lumbir of 433 people/km*. Also, the health facilities owned by Banyumas Regency are two hospitals and 
39 health centers, a Community Health Laboratory, Community Eye Health Center, Community Lung Health 
Center, as well as Kartini Mother and Child Public Health Center [16]. The data showed that the DHF cases 
in Banyumas from 2009-2018 were 3,610. The highest occurred in 2016 with 990 cases and the lowest 
occurred in 2018 with 57 cases. The highest incidence rate was in 2016 at 55/1000 population and the lowest 
was in 2009 at 2/1000. In 2018, 57 dengue cases were reported in Banyumas, and have spread throughout the 
area. The following is a description of the 57 cases DHF in 2018. The gender characteristics of the patients 
are almost the same between men and women, with a slight difference. Most of patients were women (54%) 
as shown in Table 1. 


Table 1. Characteristics of DHF patients 


Characteristics Statistics 
Age Average 25 years old 
Youngest 3 month 
The oldest 81 years old 
Characteristics Total Percentage (%) 
Male 26 46 
Female 3] 54 
Total 57 100 


DHF cases in Banyumas have an upward trend from September to December. The peak occurred in 
December with 17 cases (29.8%). The two week graph showed that there are always DHF cases reported, 
which showed that the fever occur for two weeks as shown in Figure 1. 

The following is the DHF cases distribution at Banyumas Regency in 2018. There are 57 cases 
dominated by the areas that are close to the central government in Purwokerto. Below is a map of the cases 
distribution. Figure 2 showed that the distribution of dengue cases was mostly in three sub-districts, namely 
East and West Purwokerto, as well as Somogede. This distribution area was obtained by calculating the area 
based on the outermost points of the DHF case. 
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Figure 1. 2-Weekly graph of DHF 
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Figure 2. Distribution of DHF cases 


Figure 3 shows that DHF cases were mostly found in areas with moderate rainfall. In areas with 
fewer rainfalls, the incidence tends to be less. The areas with light or fewer rainfalls include Kebasen, 
Banyumas, Somagede, Kemranjen, Rawalo, Jatilawang, Wangon, Ajibarang, Gumelar and Lumbir Sub- 
districts. In addition, the urban area in Banyumas, around Purwokerto has moderate rainfall. Buffering on this 
map is the location of health facilities, which include hospitals, health centers, and polyclinics. 


Int. J. Public Health Sci., Vol. 10, No. 2, June 2021: 231 — 240 


Int. J. Public Health Sci. ISSN: 2252-8806 O 235 


108" 66°30"E 108°Cor"eE 109" 16°30°E 109’ 270°E 


7713'30'S 


713'30'S 


A Rainfall Station 


© Cases 


[| _} Low : 5-20 
[) Middle - 20-50 
GE Heavy : 50-100 
GN Very Heavy : >100 


7°240"S 
7°240"S 


7°3430°S 
7°3430°S 


109°SO0°E 109° 16°30"E 109° 270°E 


Figure 3. Rainfall overlay with DHF cases 


The results of the analysis showed that most of the DHF cases are located close to health facilities, 
which is within a radius of less than 1 km as shown in Figure 4. The Figure 5 shows that the distribution of 
DHE cases did not approach the river flow. There are many rivers in Banyumas, and the cases do not occur in 
river basins. Population density was obtained by finding the sub-district area divided by the total population. 
The darker blue color indicates a higher density. The following is a buffering map for dengue cases with 
rivers at Banyumas. 
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Figure 4. Buffering map of DHF cases with health facilities in Banyumas Regency 
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Figure 5. Buffering Map of DHF Cases with rivers in Banyumas Regency 


Figure 6 reveals that the DHF cases occurred in the more densely populated sub-district. In sparsely 
populated areas, there are fewer cases. Meanwhile, Purwokerto urban area which includes North, South, East, 
and West has a higher population density and also more dengue cases. Areas with the highest population 
density (>1250 people/km?) are in North, South, West, and East Purwokerto, Kembaran, Sokaraja, Wangon, 
Ajibarang, Sumbang as well as Karanglewas Subdistrict. Meanwhile, the areas with the lowest population 
density (<500 people/km?) are Lumbir and Gumelar Subdistrict. Spatially, it showed that the highest cases 
occurred in urban areas, namely North, South, West, and East Purwokerto Subdistrict with high population 
density rates. In addition, the areas with the lowest cases, namely Lumbir, have low population density. 
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Figure 6. Overlay map of DHF cases with population density in Banyumas Regency in 2018 
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The following are the color symbols in the land use map in Banyumas. The colors are gray 
(buildings), green (forests or hills), brown (field land), pink (settlements), blue (swamps) and orange (rice 
fields). The map showed that the cases were mostly in residential areas. In green areas, such as hills and 
forests, 1t was reported that the cases are rarely found as shown in Figure 7. The clustering analysis results of 
DHF cases at Banyumas Regency in 2018 as shown in Figure 8. 
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Figure 7. Overlay map of DHF cases with land use at Banyumas Regency in 2018 
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Figure 8. The clustering analysis results of DHF cases at Banyumas Regency in 2018 


The analysis results with Satscan showed that the clustering of dengue cases resulted in the most 
likely cluster (primary cluster) and three secondary clusters. Most likely clusters (primary clusters) occur on 
March 13 to April 23, 2018 with six cases and centered on the coordinate point (-7.380000 S, 109.273611 E) 
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with a grouping radius of 3.53 km and statistically significant (p=0.001). This primary cluster area is a 
combination of four administrative regions, which are Kembaran, Sumbang, Baturaden and North 
Purwokerto Sub-districts. Based on the Satscan results, it can be concluded that the distribution of DHF cases 
in Banyumas Regency is not random. This explains that as long as there is one dengue patient in one area, 
residents who lives within a radius of 3.53 km have a very high chance (high risk) of being infected. 

The first secondary case grouping in East and North Purwokerto Sub-districts was centered on the 
coordinate point (-7.418889 S, 109.22944 E) with a grouping radius of 1.82 km. As many as five cases 
occurred on June 19, to August 13, 2018 and were not statistically significant p=0.200. 

The second secondary case grouping in Patikraja Subdistrict was centered on the coordinate point 
(-7.491111 S, 109.23527 E) with a grouping radius of 0.93 km. As many as two cases occurred on August 14 
to September 10, 2018 and were not statistically significant p=0.375. 

The third secondary case grouping in Karanglewas Subdistrict was centered on the coordinate 
point (-7.395556 S, 109.180000 E) with a grouping radius of 3.33 km. As many as three cases occurred on 
September 25 to November 5, 2018 and were not statistically significant p=0.540. 


4. DISCUSSION 

There was a clustering of dengue cases, which consisted of one primary and three secondary 
clusters. The primary cluster location covers 4 sub-districts, namely Kembaran, Sumbang, Baturaden, and 
North Purwokerto. This showed that there is dengue transmission in the area. Also, as long as there is one 
patient in an area, people who live within a radius of three kilometres have a risk of contracting the disease. 
The areas identified as primary clusters are a combination of urban and suburban areas. The urban areas are 
Synonymous with densely populated areas; therefore when there are people with infectious diseases, the 
possibility of spreading it will increase. Furthermore, the mobility of population from one area to another has 
further increased the spread of infectious diseases. The results of research in Rembang, Indonesia showed 
that the mobility of population in sub-district or regency is a potential for dengue transmission [17]. Research 
in Italy showed that foreign tourists traveling for a relatively long time need to carry out serological 
surveillance related to the spread of dengue. This is because there is no licensed vaccine used to prevent the 
disease [18]. One of the risk factors for secondary infection that causes DHF is a history of having been 
previously exposed and migration to urban areas [19]. Similar results in Amazonia showed that villagers 
visiting urban areas are at risk for contracting dengue [20]. Meanwhile, primary clusters in Banyumas were 
detected from March to April. In Bali, strong seasonality of dengue incidence was observed during January to 
May [21] and Southeast Sulawesi Indonesia, the highest rate of dengue cases was found in January- March 
[22]. This result is different from the study in Pakistan which showed that the primary cluster was detected 
from August to October and the secondary was detected in September [23]. 

In 2018, DHF cases spread in almost all areas of Banyumas Regency, and the cases tended to 
increase compared to 2016-2017. This Regency is one of the areas in Central Java Province which is 
generally considered a rural area. Furthermore, this region is still dominated by active agricultural areas and 
part of the community depends on agricultural sector for their livelihoods. However, there are several sub- 
districts which are urban areas, especially in Purwokerto city and its surroundings. Purwokerto is the place 
for regional government offices; therefore this area is dominated by fairly complete offices and 
infrastructure. Recent studies showed that the incidence of DHF occurred in densely populated urban and in 
rural areas. 

Overlay of DHF cases in Banyumas Regency uses land use variables, population density, and 
rainfall. The analysis results showed that the cases were mostly in relatively dense residential areas. In urban 
areas, Aedes aegypti is more common, while in rural, Aedes albopictus is more common, this is related to 
land use. The distribution of Aedes aegypti and Aedes albopictus in urban areas is found in areas with 
residential land use. Meanwhile, the distribution of Aedes aegypti and Aedes albopictus in rural areas is found 
in areas with residential land use surrounded by gardens [24]. 

The spread of DHF is mostly found in densely populated areas. Meanwhile, density is not a major 
cause of the incidence, but it is an important risk factor in its development. The high population density and 
the close distance between houses causes the spread of DHF to be more intensive in urban areas. This is 
related to the flight distance of mosquitoes up to an average of 100 meters [25]. Furthermore, the denser 
population density in an area causes transmission to be easier and faster. In fact, the higher the population 
mobility, the more the transmission from one place to another. This is because dengue virus can be 
transmitted quickly as long as the place is an endemic area [26]. Statistical test results also showed that there 
is a relationship between population density and the incidence of DHF. Areas that are densely populated tend 
to have more cases than less dense areas [27]. 
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DHE cases in Banyumas are more common in areas with moderate rainfall, while areas with little 
rainfall have few cases. This result is consistent with a study in Jakarta which showed a relationship between 
rainfall and the incidence of dengue fever. When the rainfall is high, the hourly perching rate of Aedes 
mosquitoes also increases. During the rainy season, water reservoirs increase and are used by Aedes sp. to lay 
eggs [28]. 

In Thailand, dengue outbreaks coincide with the rainy season, mean rainfall and maximum 
temperature were positively associated with dengue incidence [29]. Therefore, climate is very important in 
the transmission of vector-borne diseases such as dengue fever, especially in the highlands of tropics. 
Although rainfall has no direct effect, it plays a role in transmission. The presence of rain forms a collection 
of water which can serve as a breeding place for Aedes sp. Furthermore, high rainfall increases humidity and 
extend the life of adult mosquitoes. Rainfall can either increase or eliminate the number of breeding places, 
and also affects the temperature and relative humidity of the air. Therefore, rainfall is an important climate 
element in vector density because high levels can cause standing water which has the potential as a breeding 
site for Aedes Aegepty mosquitoes [25]. 

Buffering of DHF cases based on the existence of health services or public health centers showed 
that most cases are in areas close to health services in the form of hospitals (less than one kilometre). The 
average, minimum, and maximum distances of DHF cases from the health center are 1,613 meters, 263 
meters, and 5,625 meters, respectively. The results of previous research showed that some of the cases were 
located within four kilometres radius of Margono Soekardjo Hospital [30]. The study in Semarang stated that 
there was a relationship between access to health services and the incidence of death due to dengue. Those 
that died likely have difficult access to health services, while those that are alive have easy access to health 
services [31]. The buffering results of the cases on rivers showed that the distribution does not approach the 
river flow. The study in Rembang, Indonesia stated that DHF cases located in zone buffer 0-100 m from 
puddle [17]. 


5. CONCLUSION 

One primary and three secondary clusters of dengue were detected in Banyumas Regency in 2018. 
The primary cluster occurred in March-April, involving four sub-districts in urban areas, namely Kembaran, 
Sumbang, Baturaden, and North Purwokerto. Therefore, the significantly identified clusters indicate dengue 
transmission in Banyumas Regency with a radius of three kilometres. 
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